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Ultracold Fermi-gases near the unitary limit are studied in the framework of Ef- 
fective Field Theory. It is shown that, while one can obtain a reasonable description 
of the universal proportionality constants both in the narrow and the broad Fesh- 
bach resonance limits, the requirement of the reparametrisation invariance leads to 
appearance of the three body forces needed to cancel the otherwise arising off-shell 
uncertainties. The size of the uncertainties is estimated. 



Ultracold Fermi gases have recently attracted a lot of attention ( see [l| and references 
therein) due to exiting possibility of tuning the strength of the fermion-fermion interaction 
through the Feshbach resonances so that the scattering length may become much larger then 
the typical scale in the system. This large scattering length is the main dynamical factor 
resulting in establishing the so called unitary limit (UL) which is believed to be universal in 
a sense that the only relevant energy scale is given by that of the non-interacting Fermi gas 

3 k 2 

E GS = iE FG = i-^ (1) 

where M and kp are the fermion mass and Fermi momentum correspondingly and £ is the 
universal proportionality constant, which does not depend on the details of the interaction. 
The other dimensional characteristics of the cold Fermi-gas in the UL such as paring energy 
or chemical potential can also be represented in the same way. 

Large scattering length implies nonperturbative treatment. The most "direct" nonper- 
turbative method is based on the fixed-node MC approach . However, being potentially 
the most powerful calculational tool, direct numerical simulations still have many limitations 
related to finite size effects, discretization errors, trial wave function dependence etc which 
may even become amplified in certain physical situations (the system of several fermion 
species is one possible example). All that makes the development of the analytic approaches 
indispensable. Several such approaches have been suggested so far ranging from the exact 
renormalisation group 0, Q and expansion in terms of dimensionality of space [6| to more 
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phenomenological approaches using the density functional method [7[ and many body vari- 
ational formalism 8| . The "world average" for the value of £ is 0.42 ±0.002. It is important 
to emphasise that UL refers to the idealised situation with an infinite scattering length and 
vanishing effective radius which is the case of a broad Feshbach resonance. Even small but 
finite effective radius sets the other scale so that the system may deviate from the strict 
UL. This is the case in the narrow Feschbach resonance limit and also in nuclear /neutron 
matter where the experimental value of the effective radius is only one order of magnitude 
smaller then the scattering length and should therefore be taken into account. The purpose 
of this paper is to analyse the system of the cold Fermi atoms in the both narrow and broad 
Feschbach resonance limits in the framework of effective field theory (EFT) and study the 
general constraints EFT imposes on theoretical approaches describing the system of cold 
Fermi atoms in both UL and nearby. 

EFT is based on the fact that the low scale dynamics is only weakly dependent on 
the details of the interaction at small distances. The low scale phenomena can then be 
described by a local effective lagrangian with some effective coupling constants reflecting 
the short range dynamics in some effective, indirect way. The physical amplitudes which 
can be derived from this lagrangian take the general form of the expansion in powers of the 
low scales involved with implicit assumption than all low scale are "natural" in a sense that 
all the possible dimensionless ratios are of order unity. As we have already mentioned the 
main feature of the cold Fermi atoms in UL is the large value of the scattering length. It 
makes the use of the canonical EFT impossible as the large scattering length introduces a 
new "unnatural" scale to the problem so that the power expansion is no longer valid. To 
overcome this difficulty it was suggested to iterate the leading term of the interaction to 



all orders by so 
perturbation 0, 



ving the Lippmann - Schwinger (LS) equation and to treat the rest as a 
3- 



Being a proper field theory EFT must be regularised and renormalised. Besides, it must 
satisfy the reparametrisation invariance requirement which means that, although one can 
choose different representations for the field operators in the effective Lagrangian, the phys- 
ical amplitudes should remain the same in any representation. In formal field theory this 
statement is known as a equivalence theorem llj. In a more phenomenological language 
it means that the on-shell observables must be independent on the parametrisation used 
for the off-shell part of the fermion-fermion interaction. In the context of many fermion 
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systems with arbitrary scattering length the physical consequences of the reparametrisation 



invariance were considered in [12J. In this paper we use the findings of |12J to analyse the 
system of cold Fermi atoms at and around UL. The other important general requirement, 
usually called renormalisation group (RG) invariance, is the independence of the on-shell 
physics on the renormalisation parameters like cutoff or subtraction point. To comply with 
RG invariance the fermion-fermion scattering amplitude must satisfy the RG equation. 

According to EFT the physical amplitudes at low scale can be derived from the effective 
Lagrangian with purely short range interactions. The corresponding LS equation for the 
fully off-shell amplitude T(k', k,p) takes the form 

T(k',k,p) = V(k',k,p) + M [ ^V(k',q,p) 2 T(g ^5 . (2) 

J 2n z p z — q z + it 

Here we use k and k! to denote relative momenta and the energy dependence is given 

by p = y/ME, the on-shell momentum corresponding to the centre-of-mass energy E. One 

possible form of the interaction can be written as 

V = Vt = C + C' 2 p 2 . (3) 

Since it depends on on-shell momentum the interaction is purely energy dependent. The 
interaction can be written in a separable form so that the LS equation can easily be solved 
analytically. The resulting T-matrix takes the form 

1 (C 2 I 3 - l) 2 



T(p) C + C 2 h + k 2 C 2 (2-C 2 I 3 
where the loop integrals are 



Hp), (4) 



M 



dqq n -\ (5) 



and 



i(p) = 7& hb-r-z- (6) 



— [dq ^ 

2ir 2 J p 2 — q 2 

These loop integrals are divergent and therefore the procedure of regularisation and renor- 
malisation must be carried out. As a side remark we note that the issue of the nonpertur- 
bative renormalisation is quite a subtle problem. In contrast to the standard perturbative 
case where the usual field theoretical methods can be used to regularise the given divergent 
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graphs and then renormalize the bare coupling constants, in the nonperturbative situation 
the renormalisation of the whole integral equation must be carried out. In the case when the 
analytic solution for the scattering amplitude can be obtained, as is the case here, the renor- 
malisation of the amplitude is a rather straightforward procedure. However, if the explicit 
solution is not possible (if the long range forces are added, for example) then the special 
care is needed to perform the renormalisation in a consistent way. In this paper we follow 
the procedure used in Ref. [3] to renormalize the effective fermion-fermion amplitude in 
vacuum. We subtract the divergent integrals at some kinematical point p 2 = —fj 2 so that all 
the couplings should now depend on the subtraction point \i to ensure that the scattering 
amplitude is /i independent. The regularised fermion-fermion amplitude has the form 

T = C °+P 2C 2 m 

1 i + £(c + p 2 ^)(ip + / .)- {) 

The same expressions can be obtained in the regularisation scheme considered in |l 01 ] . 

The coupling constants can be determined from the suitable observables. For example, in 
the renormalisation scheme adopted in [10] they can be related to the low energy observables 
as 

where a and r are the scattering length and the effective radius correspondingly. Strictly 
speaking the expression for the V± is written up to next-to-leading order term in the small 
scale according to the counting scheme developed in [h]]. The LO coupling Cq scales as 
p~ x and should therefore be treated nonperturbatively. The rest can be interpreted as a 
perturbation. However, it is rather easy to solve the whole Lippmann - Schwinger equation 
and obtain the vacuum T-matrix. It is worth mentioning that in the simplest case of 
the pointlike interactions all EFT does is just complicated way of getting the well known 
phenomenological effective range expansion. The full strength of the EFT can be easily 
realised when considering more complicated situations like, for example, the interaction with 
external currents. The straightforward generalisation of the phenomenological approaches 
may lead to the conflict with gauge invariance. The EFT solves this problem in a natural 
way treating all the interaction terms on an equal footing. That's the main motivation of 
using more general approach even in the case when it can be reduced to the well known 
phenomenological approaches. 
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It is clear that the above written purely energy dependent form of interaction is not 
unique. There may exist more general form of interaction in which the energy and momen- 
tum can be treated as formally independent variables. It has the form 



V = V 2 = C + C' 2 p A + -C 2 (k 2 + k' 2 - 2p 2 



(9) 



where k and k' denote the initial and final relative momenta of the fermions and the cou- 
pling C2 describes a purely off-shell interaction. The solution of the Lippmann - Schwinger 
equation is again straightforward and we get 

1 



To 



1 + 



2(c + P zc>; 



c 2 (k 2 



k' 2 - 2p z 



-fc 2 v-k 2 )(p 2 -k' 2 )(zp +/ i) 



(10) 



where T\ is given by Eq. (j7|). 

One can see from this equation that both amplitudes coincide on-shell so that the 
reparametrisation invariance is fulfilled. 

The situation becomes much more complicated in the presence of medium. The corre- 
sponding amplitude can be calculated by solving the Feynmann-Galitskii equation 



rpr 



V + VG*G*T 



F fiF rpm 



where G F is the in-medium one-fermion propagator 



G F (k) 



0{k-p F ) 



e( PF - k) 



k — uJk + ie k — uo k — ie 
here k = (k , k) and u k = k 2 /2M. 

For the interaction V\ the solution of this equation takes the form 
1 



(11) 



(12) 



C (ji)+p*C'M 47T 



Mjm _ M(p F + Q/2) M 



4tt 2 4tt 2 
Q^A-pl+p 2 



, Pf + Q 2 + p 

Plog — — + 

p F + Q/2 - p 



(p F + Q/2) 2 +p 2 

U & __2 D2 /A ^2 



(13) 



"° P 2 F -Q 2 /4-p 2 

where pp is the Fermi momentum and the amplitude is written for the case of the non- 
zero total momentum Q of the interacting fermion pair. Having determined the in-medium 
fermion-fermion amplitude we can extract the universal constant £ by computing the energy 
density Eos of the interacting Fermi gas and using the Eq.(l). The expression for the Eqs 
is given by 

E GS = E FG + E mt , (14) 
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where 

Eint 



2 PF i-PF-Q/2 

Q 2 dQ / p 2 dpT m {p 2 , Q) + 



^ 2 Pf ljo 

QdQ / pdpT m (p 2 , Q) (pi -p 2 - Q 2 /4) 

J PF -Q/2 



(15) 



Calculation of the universal constant £ in the genuine UL does not involve energy dependent 
part of the fermion-fermion interaction so we can drop the second term in Eq.(3). We 
obtain £(UL) ~ 0.33 in the unitary limit with vanishing effective radius. Similar results 
were obtained in [fjj] using the effective range expansion for the fermion-fermion scattering 
amplitude in vacuum. In the language of EFT it corresponds to a particular choice of the 
subtraction point fi = so for the lowest-order effective coupling we get Co(/x = 0) = 47ra/M. 
We see that the effective coupling becomes arbitrary large in the UL so that it is hard to 
extract the effective couplings from observables and to justify the EFT expansion for the 
effective lagrangian in this case. As shown in [h]] each individual graph in the sum of 
the bubble diagrams for the fermion-fermion scattering amplitude goes as (4ira/M)(iap) L , 
where L is the number of loops so that each contribution in the bubble sum is bigger then 
the preceding one. It means that there is no well defined expansion parameter and it may 
result in artificial dependence on short range physics, unacceptable situation for consistent 
EFT. We emphasise, however that there is nothing wrong in using the phenomenological 
expressions like effective range expansion for the fermion-fermion amplitude to calculate 
the values for some physical quantities like £(UL). In fact, a proper use of EFT leads 
precisely to this. Therefore, there is no wonder that EFT and phenomenological description 
lead to the similar results in this (simplest) case. One only needs to keeps in mind the 
constraints and limitations of the phenomenological approach. One side remark is that 
EFT provides a natural way of incorporating gauged interactions and external currents in 
unambiguous was, the opportunity often missing in phenomenological approaches where 
the conflict with gauge invariance is rather rule then exception. It is clear that EFT and 
phenomenological approaches are very different on this level of complexity. Therefore, it 
looks more advantageous to use EFT even in the simple cases where a suitable analytic 
parametrization of the fermion-fermion amplitude like effective range expansion can formally 
provide similar answers. 

It is important to emphasise that the particle-particle (and hole-hole) summation with 
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undressed propagators represented by the FG equation provides only the simplest many- 
body approximation to the many-body problem. The further complications come from the 
self-energy and vertex corrections involving particle-hole pairs. However, they should start 
contributing at the order, prescribed by the power counting. It seems possible that for the 
dilute systems the counting should not be much different from that suggested in 1 10] for the 
fermion-fermion amplitude in vacuum since each loop involving the at least one hole line 
leads to the contribution proportional to pp which is small so that the leading contribution 
is given by the summed particle-particle ladder in accord with the power counting from 10]. 
It is then looks conceivable that self-energy and vertex corrections will contribute at higher 
order. There is however a subtlety here. Let's consider the Eq. (13) for the elementary 
fermion-fermion amplitude and take for simplicity the case with zero total momentum. The 
corresponding expression for the T matrix can be written as 

Tr = i + &\piogg%-2 PF y (16) 

In the strict UL case the term with log does not contribute and the amplitude scale as 1/pf- 
Being combined with the loop it leads to the contribution of order one. It means that particle- 
particle and hole-hole ladders as well as particle-hole rings should in principle be treated on 
the same footing. Formulation of the power counting prescription is an open problem in this 
case and the use of the either simplest approximation like FG equation or more complicated 
approaches with self-energy and vertex corrections equally requires a further justification. If 
the in-medium power counting issue is resolved then all the improvements and corrections 
can be taken into account in a systematic way. 

The situation seems to be somewhat different for the system with the finite effective range 
(Pft ~ 1) and nonnegligible shape parameter (see below). In such a system these effects 
may lead to some additional suppression for the loop diagrams involving holes so that in 
this case counting indeed is closer to the vacuum one and the corresponding self-energy 
and vertex corrections are of higher order and thus suppressed. In this case the use of the 
FG equation seems more justified. Of course,the relative contributions of the suppressed 
and unsuppressed terms depend on a concrete value assumed for the effective radius. It is 
important however to keep the values of the effective range and shape parameter within the 
limits of applicability of the effective range expansion. We stress that these arguments are 
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rather crude and a consistent power counting is yet to be formulated. For example, it is 
not quite clear what kind of diagrams should contribute at the next-to-leading order and 
at what value of pft should the counting start deviating from the vacuum one. All these 
issues need to be clarified in the future EFT-based studies of strongly interacting Fermi 
systems (dilute atomic Fermi gases, neutron/nuclear matter etc). We also note that the 
available experimental data on cold fermionic atoms correspond to the case of the broad 
Feshbach resonance with ppr << 1 so that in this context the system of cold fermionic 
atoms with a large scattering length and the finite and modest effective radius refers to some 
future experiments. On the other hand the many-fermion system of this type is realised in 
nuclear /neutron matter. 

With the effective range effects included we have obtained the value £ = 0.48 at pft ~ 1. 
The result is rather close to the "world average" and seems to bring a certain credibility 
to the approach developed here. The results of calculations in the case of narrow Feshbach 
resonance with finite effective radius are shown on Fig.l as a function of the dimensionless 
parameter ppr. One notes, that the contribution due to the pairing interactions is known 
to be fairly small [is! ]. 

As one can see from the Fig.l the corrections are quite significant already at Pf^c ~ 1 
and grow as the value of ppr increases. We found £ ~ 0.68 in the case of neutron matter with 
scattering length and effective radius being —18.5 fm and 2.7 fm correspondingly. The large 
value of the effective radius contribution suggests that one might expect the nonnegligible 
contribution form the next term of the low energy expansion of the effective fermion-fermion 
interaction which is proportional to C^p 4 . The corresponding T-matrix in a free space takes 
the form 

1 1 M C 2 2 C' 2 fC' 2 CA 4 

where we have written the T-matrix in a form, consistent with the counting rules suggested 
in [h]]. The coupling C4 entering the effective lagrangian at next-to-next leading order can 



be related to the so called shape parameter 



101 ] . Again using the neutron matter parameters 



we found approximately £c* 4 ~ 0.73. We note that the size of the correction, while being 
non- negligible, suggests that the contributions of the higher order terms can be neglected. 

All that looks reasonable but the word of caution is in order here. Let us turn to the 
more general case of the interaction V 2 with both energy and momentum dependence. The 
Feynmann-Galitskii equation can be solved in the same way and after putting the amplitude 
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FIG. 1: The universal parameter £ as the function of pi?r. 



TV 1 on-shell we obtain 



rpm 

1 2 



rpm ofrrm^ QW M 3 



(18) 



We see that two interactions which resulted in the same physical amplitudes in vacuum, 
lead to the different on-shell T-matrices in the presence of the fermion medium so that 
the reparametrisation invariance is not satisfied. Neither is satisfied the RG invariance re- 
quirement as the renormalisation performed at different subtraction points leads to different 
results for the physical observables. In other words, the physical observables still depend 
on the off-shell behaviour assumed for the fermion-fermion interaction. This is clearly the 
unsatisfactory situation which should be corrected. The general analysis of this problem was 



given in [12] so that here we give just a summary of the main points from 12j. Firstly, the 
hint on how to cancel the unphysical contributions comes from the second term in Eq.(17) 
which is proportional to the density. The same structure arises from a three-body (3B) 
contact interaction so that the 3B forces could probably be used to achieve the required 
cancellation. Secondly, we note that the contribution of the off-shell term is driven by the 
coupling constant which cannot be extracted from any physical observables which are 
defined on-shell. In the lowest order in the off-shell coupling Gi the ground state energy is 
determined by the Hugenholtz diagrams shown in Fig. 2, where the solid dot denotes an in 
medium NN vertex and thick lines are dressed fermion propagators. Hugenholtz diagrams 
are versions of Feynman diagrams which explicitly incorporate antisymmetry of the inter- 
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FIG. 2: Hugenholtz diagrams for the ground state energy at first order order in C2 (the open 
triangle). 




FIG. 3: Hugenholtz diagrams for the ground state energy at first order in the three-body force 
(the open square). 

actions. Internal lines represent Feynman propagators which describe both particles and 
holes. The arrows represent the flow of quantum numbers such as baryon number. Each 
topologically distinct diagram should be multiplied by a symmetry factor to take account 
of the number of ways it can be constructed from the antisymmetric vertices. More details 
of these diagrams and the rules for evaluating them can be found in the textbooks [if], fl7 |- 
The diagrams in Fig. [2] give rise to many different contributions, which can be identified 
by iterating the equations for the in-medium NN vertex and dressed propagator. The ap- 
proximations commonly used in many-body physics, typically amount to replacing the full 
NN vertex by a G- or T-matrix including both particle-particle {p — p) and hole- hole (h — h) 
adders. Including ph rings as well as pp and hh ladders leads to the parquet approximation 



16]. 



As shown in |12| the terms containing the unwanted off-shell contributions can indeed be 
exactly canceled against the contributions of a contact 3B interactions with three distinct 
topological structures shown in Fig. 3. In the simplest Brueckner-Hartree-Fock approxima- 
tion {la . \v\ . in which propagators are dressed and the in-medium NN vertex is obtained 
by iterating the potential in the p — p and h — h channels the unphysical contributions can 
be shown to cancel with Fig. 3 (a). When p — h channel is added so that one iterates the 
interaction in all {p — p, h — h and p — h) channels all three graphs from Fig. (3) are needed 
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to achieve the required cancellations. 

We note that it is rather hard to make a rigorous statement on whether it is at all 
possible to formulate the general approach based on the 2B forces without the 3B ones and 
simultaneously satisfying the reparametrisation invariance requirement. The form of the 
second term in Eq.(17) seems to indicate that this is not possible but this is admittedly 
suggestive argument. In the context of the above discussion the more rigorous statement 
would be that for the most popular and widely used many-body approaches such as p—pj h—h 



ladder, parquet and even advanced parquet [16| approximations inclusion of the 3B forces is 



required to satisfy the reparametrisation invariance theorem. An additional support for this 
statement comes from the EFT studies of the few-body systems [3] where the 3B forces 
are needed to carry out a consistent renormalisation procedure. Moreover, as shown in [3] 
the corresponding 3B vertex must be promoted to a leading order in effective Lagrangian. 
All that strongly suggests that the 3B forces must necessarily be included at any level of 
truncations used so far in theoretical calculations. 

The importance of the higher order diagrams with both the 2B and 3B forces could be 
estimated more quantitatively from some power counting rules. Unfortunately, as we pointed 
out above, establishing such a counting for the strongly interacting Fermi system is still an 
open and very challenging problem. 

As we already mentioned, the off-shell parameters cannot be extracted from the on-shell 
physics so one can only theoretically estimate the off-shell contribution and, hence, the 
strength of the 3B forces needed to cancel it. The possible estimate could be based on the 
assumption that the term with the coupling Ci gives the contributions of the same order as 
those related to its on-shell "cousin" C' 2 . It leads to 

C 2 ~ C' 2 ~ 0(1), (19) 

if n ~ pp. One notes that the 3B forces will also depend on the subtraction point // to satisfy 
the RG invariance. Of course, the estimates obtained for the 3B forces are very crude and 
much quantitative work remains to be done to properly take their effect into account. Apart 
from the "destructive" role of cancelling the off-shell contributions the 3B forces should also 
play a "constructive" role in bringing the theory to better agreement with the data. The 
examples include the binding energy of the triton or low energy neutron-deuteron scattering 
[3]. The other example, which is very general and hold for any three body system with 
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infinite scattering length in the two body subsystems is the Efimov effect [^(J which states 
that in such a system there exist infinitely many three-body bound states. It clearly should 
have a huge influence on the dynamics of the cold Fermi gases near the unitary limit and 
should therefore be taken into account in any consistent theoretical approach. Work in this 
direction is in progress. 

The bottom line of this discussion is that, in spite of the fact that the EFT motivated 
studies of the Fermi gases near the unitary limit may result in the reasonable numbers for 
the physical observables they should be interpreted with a caution as neither power counting 
nor reparametrisation (and renormalisation) invariance issue is satisfactory implemented in 
theoretical schemes at present. 



[1] E. Timmermans et al., Phys. Lett. A 285, 228 (2001). 

[2] J. Carlson et al., Phys. Rev. Lett. 91, 050401 (2003); 

[3] G. E. Astrakharchik et al., Phys. Rev. Lett 93, 200404 (2004). 

[4] B. Krippa. larXiv:0704.3984l 

[5] S. Diehl, H. Gies, J. M. Pavlovski and C. Wetterich, [cond-mat/0701198 



[6] Y.Nishida and D. T. Son, Phys. Rev. Lett. 97, 050403 (2006); Y.Nishida and D. T. Son, 



cond-ma t/0607835; T. Schafer, cond-m at/0701251 P. Arnold, J. E. Drut and D. T. Son, 



cond-mat/0608477 



[7] A. Bulgac, cond-mat/0703526 



[8] R. Haussmann et al., Phys. Rev. A 75, 023610 (2007). 

[9] S. Weinberg, Nucl. Phys. B363, 3 (1991). 

[10] D. B. Kaplan, M. J. Savage and M. B. Wise, Nucl. Phys. B534, 329 (1998). 

[11] R. Haag, Phys. Rev. B112, 669 (1958). 

[12] M. C. Birse, B. Krippa, N. R. Walet and J. A. McGovern, Phys.Rev.C67, 031301 (2003). 

[13] J. Gegelia, |nucl-th/9802"038| 

[14] H. Heiselberg, Phys.Rev.A63, 043606 (2001). 

[15] O. Elgar0y and M. Hjorth-Jensen, Phys. Rev. C57, 1174 (1998). 

[16] J. P. Blaizot and G. Ripka, Quantum theory of finite systems. (MIT Press, Cambridge, 1985). 

[17] J. W. Negele and H. Orland, Quantum many-particle systems. (Perseus, New York, 1988). 



13 



[18] P. F. Bedaque and U.van Kolck, |nucl-th /0203055 , Ann.Rev.Nucl.Part.Sci.52 339, (2002). 
[19] A. Stadler et al, Phys.Rev.C51, 2896 (1995). 
[20] V. N. Efimov, Phys. Lett. 3B, 563 (1970). 



